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Abstract. It has been reported that traveling waves propagate periodically and stably in sub-excitable 
systems driven by noise [Phys. Rev. Lett. 88, 138301 (2002)]. As a further investigation, here we observe 
different types of traveling waves under different noises and periodic forces, using a simplified Oregonator 
model. Depending on different noises and periodic forces, we have observed different types of wave prop- 
agation (or their disappearance). Moreover, the reversal phenomena are observed in this system based on 
the numerical experiments in the one-dimensional space. As an explanation, we regard it as the effect of 
periodic forces. Thus, we give qualitative explanations to how reversal phenomena stably appear, which 
seem to arise from the mixing function of the periodic force and the noise. And the output period and 
three velocities (the normal, the positive and the negative) of the travelling waves are defined and their 
relationship with the periodic forces, along with the types of waves, are also studied in sub-excitable system 
under a fixed noise intensity. 
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1 Introduction 

The effects of noise on nonlinear systems are the subject of 
intense experimental and theoretical investigations. Noise 
can induce transition [112] , bifurcations [3] , and stochastic 
resonance |4|5|6|7j . Especially, in Ref. [7] the synchroniza- 
tion of spatiotemporal patterns were observed in an ex- 
citable medium by the numerical evidence. Moreover noise 
can enhance propagation in arrays of coupled bistable 
oscillators {6)9110)11] ■ In an excitable system, an exter- 
nal periodic forcing can dramatically change its behav- 
ior. As reported by previous documents, the phase lock- 
ing, quasi-periodicity, period doubling, and chaos were ob- 
served [12]. The temporal evolution of the concentration 
patterns has been modeled by partial differential reaction- 
diffusion equations. Such models include oscillatory, ex- 
citable or bistable systems with either none, one or two 
linearly stable homogeneous states [13114] , It is also well 
known that in sub-excitable systems noise also can in- 
duce travelling waves [15 , drive avalanche behavior [16] . 
and sustain pulsating patterns and global oscillations |17j . 
Sub-excitable systems under noises and periodic forcing 
are able to send out travelling and spiral waves. Belousov- 
Zhabotinsky (BZ) reaction [18119] is a popular symbol in 
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nonlinear dynamical realm to study excitable and sub- 
excitable system. It has been widely agreed that the noise 
and periodic forcing play a very important role in wave 
propagation and stability. 

Recently, it was observed that noise can support wave 
propagations in sub-excitable [20117115] systems due to 
a noise-induced transition [21122] . In these studied, the 
medias are static, and transports are governed by diffu- 
sion [23]. However in many situations, the medias are not 
static but subject to a motion. For example, stirred by a 
flow, or by the periodic forcing, the convective-like phe- 
nomena were observed due to electric field in Ref. [24], 
which occurs especially in chemical reactions in a fluid en- 
vironment. In such cases, diffusive transport usually dom- 
inates only at small spatial scales while mixing due to the 
flow in much faster at large scale. In Refs. [25126127] , the 
authors shown that in an inhomogeneous self-sustained 
oscillatory media, an increasing rate of mixing can lead to 
a transition to a global synchronization of the whole me- 
dia. Especially, Ref. [28] showed that the interplay among 
excitability, noise, diffusion and mixing can generate vari- 
ous pattern formation in a 2D FigzHugh-Nagumo (FHN) 
model subject to the advection by a chaotic flow. Here, 
we research the effect of noise and periodic forcing on 
sub-excitable systems using the Oregonator model in one 
dimension, which advances from the BZ reaction. The re- 
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versal phenomenon is not observed in the previous doc- 
uments about the propagation of travelling waves. It is 
founded in this paper, which relates to a new concept. In 
Ref. 30J, Richard A. Albanese refers to reversal concept in 
wave propagation concern extraction of information about 
distant structural features from the measurements of scat- 
tered waves, but it is irrelevant the excitable system. 

In our paper, we define that the general traveling waves 
propagate forward in one constant direction and vice versa. 
Under this definition, we found the reversal phenomena in 
our numerical simulations. However in our simulation we 
have discovered that after some time, the waves change 
its propagation direction and turn backward to travel in 
the opposite direction. The traveling waves propagate for- 
ward and backward alternately and periodically. That is 
called the reversal phenomenon (see the supplementary 
material on-line movie for this phenomenon, Movie- 0). In 
the mathematical language, the definition of definition of 
an "reversal phenomenon" as at time t, the spatial posi- 
tion of a traveling wave front is at point x. After some 
time At > 0, the wave front reaches point x again with 
reversal direction. In fact, this phenomenon was observed 
in the FHN model by numerical simulations [28]. In this 
article, our focuses are the effect of the periodic forcing 
and noise on the propagation of traveling waves in sub- 
excitable systems. 



2 Model 

Most of the systems we are interested in reside in a d- 
dimensional world. This means that our variables (fields 
or concentrations) depend on time and space. In present 
paper, the starting deterministic model [5D] is based on 
partial differential equations, and when the randomness 
is introduced we transform them into stochastic partial 
differential equation. A representative example is the de- 
terministic reaction-diffusion equation, 



^(x,t) 
dt 



f(4>(^t)^)+vv 2 4>(^t), 



(i) 



where </>(x, t) represent the density of a physical observ- 
able, /(0(x, i),/i) is a nonlinear function of the field <f> 
and ji denotes the relevant control parameter. The above 
equation can be made more complicated when considering 
vector fields, higher-order derivatives, or nonlocal opera- 
tors. The effect of fluctuations is introduced through a 
stochastic process or noise £(x, t) with well controlled sta- 
tistical properties. As a result, we expect that the new 
equation governing our system will have the generic form 



d<t>(*,t) 
dt 



/(</>(x, t), n) + 2?V 2 </>(x, t) + g(4>)£(x, t). (2) 



We take into account this standard example of stochas- 
tic partial differential equation and the two- variable Oreg- 
onator model [19135] that is famous for its convenience 
to study the property of diffusion-reaction systems. Our 



modified model adds both noise and the periodic forces, 
which is, 



^ = -f(u,v) + D u V 2 u+(D s at) 
h £ 

-^=g(u,v) + D v V 2 v 



(3a) 

dx 

(D s at)+E(t))£, (3b) 



where f(u, v) — u(l — u) — fv-j^L, and g(u, v) — u — v. V 2 
is the Laplacian operator in Cartesian coordinates, u and 
v represent the concentration of HBr02 and the catalyst 
2Ce + , respectively. Here, the external electric field can 
be considered as a spatially uniform electric field, which 
includes two parts: stochastic forcing, £(t) (we use the 
notation of for irrelevant in space) and the periodic 
force , E(t). Both of them depend on time t with the form 
of Gaussian noise and with the periodic function respec- 
tively, denotes Gaussian white noise with = 
and (£(t)£(i )) = 2D s S(t — t'), a typical temporally varied 
Gaussian white noise and here D s is the intensity of noise. 
E(t) is the periodic force, where the sine periodic force 
is chosen, E(t) = Fsm(^-t). F and T- ln are the intensity 
of the periodic forcing and the input period, respectively. 
The effect of electric field is convective-like, as discussed 
in Ref. [24]. D u and D v are the dimensionless diffusion 
coefficients of u and v, D u — 1.0, D v — 0.6 [3"T] . 

The dynamical system ([3]) is simulated in one-dimensional 
space by Euler-Maruyama Method [32] with zero flux bound- 
ary conditions on a space length of 500 elements. The 
space step is Ax = 0.15 space unit and the time step 
is At = 10 -3 time unit. These parameters are chosen to 
make the simulation process relatively stable and the in- 
formation of wave propagation can be relatively all-around. 
In order to avoid the simulation going to the negative re- 
gion of u, we let u > q, if u < q, u — q, according to the 
method of the Refs. [33134] . The simulations have been 
done as follows: we excite 3 elements at the left bound- 
ary in the our systems under sub-excitable state in the 
one-dimensional space, which serves as the wave source. 
The leapfrog method for the advection terms; an implicit 
method for the diffusion terms; and a simple explicit Euler 
method for the reaction terms. Letting [7™ ~ u(xj,t n ) and 



V n w v(xj, t n ), where x 



i+i' 



= Ax and t n+1 -t n = At, 



results in the following discretised system: 



U, 



n+l 



UV 



At 



yn+1 _ yn 
3 3 



-f(U?,V?) + D u 



U?+? - 2U} 1+1 + Upi 



n+l 



rn+1. 



At 



+(D£(t) + E(t)) 
g(JJ?,V?)+D v 
+ (D s £(t) + E(t)) 



(Ax) 

r^+i-tT-i 



2Aa 



(4) 



Vj$ L 1 -2Vp +1 + V? +l 



(Ax) 2 

rV n — V n 



2Ax 



(5) 



where D s , F, and T- ln are the control parameters. Note that 
here the noise term £(t) is Gaussian and white. This is a 
very reasonable assumption for internal noise, which rep- 
resents many irrelevant degrees of freedom evolving very 
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Fig. 1. (Color online) The phase diagram of the e versus /. 



short temporal and spatial scales. It has a probability den- 
sity function of the normal distribution (also known as 
Gaussian distribution). In other words, the values that 
the noise can take on are Gaussian distributed. 



3 Result 

Extensive testing is performed through numerical simula- 
tions to the described model ([3]) and the qualitative results 
are shown this section. 



3.1 Sub-excitation and reversal phenomena 

q, e and / are parameters related to the BZ kinetics, de- 
termining the sub-excitation of the system, q — 0.002 (35] . 
First we simulate the dynamical system (|3]) without noise 
and periodic force to confirm the sub-excitable region. The 
result is shown in Fig. Q] In region (A) no traveling waves 
are produced in the system; in region (B) traveling waves 
are sent out and propagating in the system; region (C) is 
not included in the excitable confine. In the area between 
region (A) and (B), traveling waves are sent out but die 
away gradually, which is so-called the sub-excitable re- 
gion, corresponding to the blue line in Fig. [T] Here we 
set e = 0.1 and / = 2.435 to fix our system into in a 
sub-excitable region. When the parameters e and / are 
deeply in the sub-excitable domain, the same results are 
observed by the numerical experiments, such as, e = 0.025, 
/ = 2.435 and e = 0.025, / = 4.0. 

It turns out that the influence of noise is rather impor- 
tant in the sub-excitable media. There is an extreme case 
for the system ((3|), that is, when the intensity of noise 
equals to 0, i.e. D s = 0, and only the periodic force is 
present, the different types of wave propagation are sum- 
marized in the Table [TJ Table [T] reveals that there exists 
critical input period T* n for each intensity of periodic force 
F (F > 10.0), so that the system creates traveling waves 



Fig. 2. (Color online) The spatiotemporal plot of the vari- 
able u for the system ((3j), where the ordinate is time evolu- 
tion and the abscissa represents the spatial location. The white 
part indicates wave front. (A) only noise present, D s — 14.0, 
and F = (see the Movie-2, additional movies format avail- 
able from Journal Web); (B) Both noise and periodic forcing 
present, D s = 10.0, F = 10.0, and T in = 8.0 (see the Movie- 
3, additional movies format available from Journal Web); (C) 
Both noise and periodic forcing present, D s — 14.0, F = 10.0, 
and Ti n = 8.0 (see the Movie-4, additional movies format avail- 
able from Journal Web). 



periodically and waves propagate persistently. If F < 10.0, 
the system can only send out one wave, and there is no 
critical input period T* n ; The system sends out traveling 
waves periodically but the waves disappear quickly when 
the Ti n > T* n (see Movie-1 for this case, where D s = 0, 
T in = 8.0, and F = 10.0); when F = 10.0, T* n = 5.0; 
F = 20.0, T* n = 4.0, 5.0; F = 30.0, T* n = 3.0, 4.0, 5.0. 

Now, we turn on the noise and the periodic forces in 
the system ([3]) and investigate their effects when the in- 
tensity of noise takes different values. Fig. [2] shows that 
traveling waves propagate with different noise and peri- 
odic forces. The abscissa is the spatial location and the 
ordinate is the evolution of time. The white part indi- 
cates wave crests. As shown in Fig. EJA), if only noise 
present, traveling waves are produced irregularly but they 
can propagate stably [see Movie-2]. If periodic force and 
noise both present, travelling waves are produced period- 
ically and reversal phenomena appear, however traveling 
waves die out quickly when the intensity of noise is small 
[see Movie-3]. The travelling waves are produced periodi- 
cally and they propagate stably, and reversal phenomena 
also appear when the intensity of noise is increased [see 
Movie-4), as shown in Fig. [2IC). In Fig. [3] we show the 
phase diagram about the reversal phenomena with respect 
to the D s -F parameter space, in which the reversal trav- 
elling waves will die out quickly within the region I but 
propagate persistently within the region II. One can see 
that the reversal waves sensitively depend on the intensity 
of periodic forces and the noise intensities for the fixing 
Ti n . For example (see Fig. [3|), the emergence of reversal 
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Table 1. Traveling waves propagation under different F and T n with noise intensity D s — 0; symbol —1 indicates the system 
sends out only one wave; symbol — indicates the system sends out traveling waves periodically but the waves disappear quickly; 
symbol + indicates the system sends out traveling waves periodically and the waves can propagate persistently. 





Tin < 2.0 


T n = 2.0 


T„ = 3.0 Tin = 4.3 


Tin = 5.0 Tin = 6.0 Tin > 6.0 
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-1 
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+ - - 
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-1 


-1 
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Fig. 3. (Color online) Phase diagram about the reversal phe- 
nomena with respect to D a -F parameter space. The parame- 
ters are the same as Fig. [2] but T n = 8.0. The reversal waves 
will die out quickly within region I and propagate persistently 
within region II. Note that the log scale for T n . (T, the critical 
values from simulation results.) 



waves is a sensitive relationship at before and after the 
certain critical value. For the large F, the curve is sharp 
decreasing, otherwise the cure is slow decreasing. So we 
can obtain conclusion from Figs. [2] and [3l that the stabi- 
lization of traveling waves are due to noise, but the pe- 
riodicity owe to the function of periodic force. Periodical 
and stable traveling waves, as well as reversal phenomena, 
are produced in sub-excitable system driven by noise and 
periodic force. 

To characterize the relation between the periodic force 
and reversal phenomena, we give some qualitative expla- 
nations. Fig. 0] shows the relationship between the prop- 
agation of traveling waves with the periodic force E(t) — 
Fsin(^r L i). The left window of Fig. [4]is the periodic force 
whose abscissa is E(t) and the right window represents 
the spatiotemporal plot of variable u whose abscissa is 
the spatial location. The left and right window share the 
same ordinate which is time. From Fig. [5] we can easily 
obtain that corresponding to each input period the sys- 
tem sends out a traveling wave, which is 1 : 1 frequency 
locking. The dotted line in the left window of Fig. 2] divide 
out the positive and negative parts of the periodic force. 




Space (a.u.) 



Fig. 4. (Color online) The propagation of traveling waves with 
the periodic force E(t) = F sm(j^-t). The ordinate is time 
evolution. For the left window the abscissa represents E(t) and 
for the right the abscissa is space location. The left and right 
window share the same ordinate which is time. The parameters 
are D s = 14.0, F = 10.0, and T in = 8.0. 




Space (a.u.) 



Fig. 5. (Color online) The propagation of travelling waves 
with the periodic forcing E(t) = E-(t) (or E(t) = E+(t)). The 
ordinate is time evolution. For the left window the abscissa 
represents E(t) and for the right the abscissa is space location. 
The left and right window share the same ordinate which is 
time. The parameters are D s = 14.0, F = 10.0, and T in = 8.0. 
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We find out that if the periodic force is negative, traveling 
waves propagate forward, whereas if the periodic force is 
positive, traveling waves propagate backward in the oppo- 
site direction. Thus reversal phenomena appear. 

For further investigation, we test several different kinds 

I 2t 

of periodic force for system (J3J). They are E n (t) = (— l) LT in- 
(the rectangle periodic force, here the [n\ denotes the 
integer of the n), E-{t) = -|Fsin(|^i) and E + (t) = 

|Fsin(|r-t)|. When the periodic forcing is E n (t), trav- 
elling waves are produced periodically and they propa- 
gate stably, and reversal phenomena appear (resemble the 
Fig. G3TJ)). When the periodic force is E_(t) or E+(t), 
there are only foundations formed but no travelling waves 
are sent out, as shown in Fig. [5] So the reversal phe- 
nomenon is due to the alternation of the positive and neg- 
ative values of the periodic force. If each value of the pe- 
riodic force was positive (or negative) , no travelling waves 
will be sent out and of course no reversal phenomena will 
appear. 



3.2 The output period and the three velocities 

In this section we focus on the effect of the periodic force 
E{t) =Fsm{^-t) on the propagation of travelling waves. 
Here, the noise intensity D s is fixed at 14.0. Furthermore, 
we define several quantities of the traveling waves. The 
output period T ou t is defined as follows: T is the time in- 
terval between the ith wave and the i + 1th wave, m waves 
are taken into account and the average value of them is 
£ T 

Tout, where Tout = ■ The three velocities (the nor- 
mal, the positive and the negative) of the traveling waves 
are defined as follows: when traveling waves propagate for- 
ward, there is a mean propagation velocity which is the 
positive velocity, and denoted by V+. The mean velocity 
of the backward propagating waves is likewise denoted by 
V-. In addition, V is defined here as the average velocity 
of the whole propagation of traveling wave. The output 
period Tout, the normal velocity V, the positive velocity 
V+ and the negative velocity V_ are shown in Fig. [6] The 
slopes of the red, blue and yellow lines in Fig. [6] are 1/V, 
1/V+ and 1/V_, respectively. 

Based on the definition, we now focus our simulations 
on the relationship between periodic force and the prop- 
erties of travelling waves. The results are shown in Fig. [JJ 
Above all we interpret how the periodic force E(t) affects 
the output period T ont . First, we fix F = 10.0 (the in- 
tensity of the periodic force) and study how T out changes 
with Ti n (the input period), as shown in Fig. [Jja). One ob- 
serves that the output period increases linearly with the 
input period. Moreover, the output period is the same as 
the input period (the slope equals to 1). Second, we fix 
T in = 8.0 and study the relation between T out and F, as 
shown in Fig. UJ^b) , from which one observes that the out- 
put period is a constant and independent of the intensity 
of the period forcing. So we can draw a conclusion that 
the output period is the same as the input period and is 




Fig. 6. (Color online) The sketch map of T out , V, V+, and V- 
with D B = 14.0, F = 10.0, T in = 8.0. 
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Fig. 7. (Color online) The effect of periodic force on the prop- 
erty of traveling waves, D s = 14.0. (a) The plot of the output 
period with respect to the input period, F = 10.0. (b) The plot 
of the output period with respect to the intensity of periodic 
force, T n = 8.0. (c) The plot of the three velocities of the trav- 
eling waves with respect to the input period, and F = 10.0. 
(d) The plot of the three velocities of the traveling waves with 
respect to the intensity of periodic force, T n = 8.0. 



independent of the intensity of the periodic force, under a 
fixed noise intensity. 

Next we investigate the influence of the periodic force 
E(t) on the three velocities of the travelling waves (the 
normal V, the positive V+ and the negative VL). First 
we fix the intensity of the period forcing F = 10.0 and 
see the changes of V, V+, and V_ with the input period 
Tn, respectively. The results are shown in Fig.UJc), from 
which we observe that the normal, positive and negative 
velocities approximate constant values, respectively, inde- 
pendent of the input period. Then we fix T; n = 8.0 and 
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study the changes of V, V+ and V- with F (see Fig.JTJd)), 
which shows that the normal velocity is constant indepen- 
dent of the intensity of the periodic forcing; the positive 
velocity increases with the intensity of periodic force, while 
the negative velocity decreases with the intensity of peri- 
odic force. So we can conclude that the normal velocity 
is independent of the periodic force; the positive the neg- 
ative velocities have the same trend with the intensity of 
periodic forces and do not depend on the input period. 

It is natural to presume from the above results that 
the three velocities are interrelated. Through quantities 
of data, we obtain the relationship among the normal, 
positive and negative velocities, as follows 

2V = V+ + V-, (6) 

that is, 

V* = V+-V = V -V-. (7) 

So the positive velocity is the normal velocity plus the 
V* , and the negative velocity is the normal velocity minus 
the V*. 

4 conclusion and discussion 

In conclusion, noise and periodic force play a very impor- 
tant role on the production, propagation, and stability of 
the traveling waves in sub-excitable systems. Noise can 
induce traveling waves to propagate stably. It can also 
support wave propagation in sub-excitable [20 17 15j me- 
dia due to a noise-induced transition [21122136] . In these 
studies, the media are static, and transport is governed 
by diffusion. In many systems, the media are not static, 
but subject to a motion, for example, when stirred by a 
flow, or by the oscillatory electric field (with period) . This 
occurs especially in chemical reactions in a fluid environ- 
ment. In this case, usually diffusive transport dominates 
only at small spatial scales while mixing due to the flow in 
much faster at large scale. In this paper, we investigate the 
sub-excitable system using a simplified Oregonator model, 
and the propagation of traveling waves in the presence of 
both noise and periodic force. Depending on noise and the 
periodic forcing we have observed different types of wave 
propagation (or its disappearance). Moreover, the rever- 
sal phenomena is observed in this system based on the 
numerical experiments in one-dimensional space, we give 
qualitative explanations to how reversal phenomena ap- 
pear, which turns out to be the periodic force. The output 
period and three velocities (the normal, the positive and 
the negative) of the traveling waves are defined and their 
relationship with the periodic forcing are also studied in 
sub-excitable system with a fixed intensity of noise. 

The periodic force can make traveling waves produced 
periodically and reversal phenomena appear. The reversal 
phenomenon results from the alternation of the positive 
and negative values of the periodic forcing and noise. For 
the special case, we give the phase diagram about the re- 
versal phenomena with respect to D S ~F parameter space, 
from which one can see that the reversal waves sensitively 



depend on the intensity of periodic forces and the noise in- 
tensities for the fixing Tj n . At last we research on the effect 
of periodic force on the sub-excitable system with a fixed 
noise intensity D s = 14.0. Under thus a condition, the in- 
fluence of periodic forces on the three different kinds of 
velocities are also acquired. The relation among the three 
velocities is 2V = V+ + V-. In Ref. [20], the authors stud- 
ied sub-excitable medium of Belousov-Zhabotinsky (BZ) 
reaction subjected to Gaussian white noise in experiments. 
They observed that at an optimal level of noise the wave 
sources of excited traveling waves become synchronous, as 
though there exists a long distance spatial correlation. 
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